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Abstract 


The origin of Jupiter and the Galilean satellite system is examined in the light of 
the new data that has been obtained by the NASA Galileo Project. In particular, special 
attention is given to a theory of satellite origin which was put forward at the start of the 
Galileo Mission and on the basis of which several predictions have now been proven 
successful (Prentice, 1996a, b, c). These predictions concern the chemical composition of 
Jupiter’s atmosphere and the physical structure of the satellites. According to the 
proposed theory of satellite origin, each of the Galilean satellites formed by chemical 
condensation and gravitational accumulation of solid grains within a concentric family of 
orbiting gas rings. These rings were cast off equatorially by the rotating proto- Jovian 
cloud (PJC) which contracted gravitationally to form Jupiter some 4 !4 billion years ago. 
The PJC formed from the gas and grains left over from the gas ring that had been shed at 
Jupiter’s orbit by the contracting proto-solar cloud (PSC). Supersonic turbulent 
convection provides the means for shedding discrete gas rings. 

The temperatures T n of the system of gas rings shed by the PSC and PJC vary 
with their respective mean orbital radii R n {n = 0,1,2,...) according as T n x R n 0 9 . If the 
planet Mercury condenses at 1 640 K, so accounting for the high density of that planet via 
a process of chemical fractionation between iron and silicates, then T n at Jupiter’s orbit is 
158 K. Only 35% of the water vapour condenses out. Thus fractionation between rock 
and ice, together with an enhancement in the abundance of solids relative to gas which 
takes place through gravitational sedimentation of solids onto the mean orbit of the gas 
ring, ensures nearly equal proportions of rock and ice in each of Ganymede and Callisto. 
Io and Europa condense above the FGO ice point and consist solely of hydrated rock 
(h-rock). The Ganymedan condensate consists of h-rock and H 2 O ice. For Callisto, NH 3 
ice makes up -5% of the condensate mass next to h-rock ( — 50%) and H 2 O ice (-45%). 

Detailed thermal and structural models for each of Europa, Ganymede and 
Callisto are constructed on the basis of the above initial bulk chemical compositions. For 
Europa (E), a predicted 2-zone model consisting of a dehydrated rock core of mass 0.912 
A/g and a 150 km thick frozen mantle of salty H 2 O yields a moment-of-inertia coefficient 
which matches the Galileo Orbiter gravity measurement. For Ganymede (G), a 3-zone 
model possessing an inner core of solid FeS and mass —0.1 16 Mq, and an outer H 2 O ice 
mantle of mass -0.502 Mq is needed to explain the gravity data. Ganymede’s native 
magnetic field was formed by thermoremanent magnetization of FejC^. A new Callisto 
(C) model is proposed consisting of a core of mass 0.826 Mq containing a uniform 
mixture of h-rock (60% by mass) and H 2 O and NH 3 ices, and capped by a mantle of pure 
ice. This model may have the capacity to yield a thin layer of liquid NH 3 - 2 H 2 O at the 
core boundary, in line with Galileo’s discovery of an induced magnetic field. 
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1. Introduction 


The successful completion of the Galileo spacecraft primary mission to Jupiter in 
December 1997 has brought a wealth of new information about this planet, its 
magnetosphere, and its family of four large Galilean satellites. While the Galileo data has 
consolidated some of our previous ideas about the makeup of the Jovian system, especially 
those which have been acquired since the Voyager 1 and 2 missions in 1979, there have 
also been many new discoveries that were completely unforeseen. For example, the find- 
ing by the Galileo Probe (Niemann et al., 1996) that the Jovian atmosphere is very dry, at 
least in the region where the Probe entered, was not predicted by any theory. Nor was the 
discovery of Ganymede’s native magnetic field, whose strength is an order of magnitude 
larger than the local ambient field (Kivelson et al., 1996). The fact that this moon is so 
strongly differentiated in its internal structure compared to its nearby twin Callisto 
(Anderson et al., 1996a, 1997b, 1998b) presents us with a deepening theoretical mystery 
(Mueller and McKinnon, 1988; McKinnon, 1997). Callisto has no native magnetic field 
and appears thermally dormant. Why are Ganymede and Callisto so different? 

With regard to Io, it now seems certain that this moon contains a partly molten 
interior (Anderson et al., 1996b), as first predicted by Peale et al. (1979). But the extent of 
the melted region and its precise chemical composition are unknown. FeS is most likely 
the chief melt constituent (Prentice and ter Haar, 1979; Prentice, 1996b). And lastly, the 
discovery of the disrupted ice crust of Europa means that throughout some portion of this 
moon’s past, it has been surrounded by a warm ocean some 100-200 km deep (Anderson et 
al., 1997a; Carr et al., 1998). Such an ocean may have formed through dehydration of an 
initial, heavily hydrated, rocky interior (Prentice and ter Haar, 1979; Ransford et al., 1981). 

This paper explores the origin, bulk chemical composition, and physical structure 
of the Galilean satellite system of Jupiter in the light of the Galileo spacecraft data. It is the 
sequel to a pre-Galileo analysis of Jupiter’s formation (Prentice, 1996a). Several important 
predictions which were put forward then have now been tested by the Galileo data. Both 
then and now it is assumed that the Galilean satellite system is a miniature version of the 
sun’s planetary system and that both systems were formed by similar mechanisms. In 
particular, the Galilean satellites, like solid planetary cores, formed through the chemical 
condensation and gravitational accumulation of myriads of grains within a concentric 
family of orbiting gas rings. These had been cast off close to the present mean orbital 


1 



positions by the primitive cloud of gas and dust which contracted gravitationally to form 
Jupiter. This modem version of Laplace’s nebular hypothesis of solar system origin is 
discussed in Section 2 below. 

In Sections 3 and 4, the gravitational contraction of the proto-solar cloud is 
examined in order to estimate the physical and chemical properties of the gas ring from 
which Jupiter’s embryonic planetary core and primitive atmosphere formed. The 
gravitational contraction of the proto-Jovian cloud is explored in Section 5, enabling a 
determination of the properties of the gas rings from which the Galilean satellites 
condensed. The initial chemical composition of each of the Galilean satellites is calculated 
in Section 6. These compositions form the basis for the construction of detailed thermal 
and structural models of the satellites. Section 7 describes the computational codes which 
generate the satellite models. Section 8 presents the results of the numerical calculations 
along with a careful comparison of these results against the Galileo Orbiter observations. 
The conclusions of this study are given in Section 9. 


2. Main Features of the Modern Laplacian Theory 

According to the Laplace (1796) nebula hypothesis of solar system origin, which is 
depicted in Figure 1, the planets accreted close to their present orbital positions from rings 
of gas that had been shed by the primitive contracting solar cloud. These rings were cast 
off by this cloud from its equatorial zone as a means for disposing of excess spin angular 
momentum during its inward radial contraction, starting from the orbit of Neptune. Of 
course, Laplace did not specify the physical mechanism for shedding discrete gas rings and 
the hypothesis was eventually abandoned since, more seriously, it could not reconcile the 
very small mass of the planetary system relative to the sun’s mass. That is, the process of 
ridding angular momentum must also lead to a substantive loss of mass from the cloud. 
This lost mass should now be shared amongst the planets. 
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Figure 1. An artist’s visualization of Laplace’s nebular hypothesis of solar system origin. 

A contracting and rotating gas cloud, originally larger in size than Neptune’s orbit, sheds 
a concentric system of orbiting gas rings. The planets later condensed somehow from 
the rings. Original drawings due to Scriven Bolton, F.R.A.S. (Whipple, 1968, Fig. 172). 

In 1928, Jeans noted that if the proto-solar cloud was very centrally condensed, i.e., 
possessed a very small axial moment-of-inertia factor / « 1, and if it rotated uniformly 
like a solid body then it could dispose of all of its excess spin momentum at the expense of 
losing little mass. That is, as Schatzmann (1949) demonstrated, if L(R e ) and M (R e ) 
denote the residual angular momentum and mass of the cloud at equatorial radius R e , then 
these quantities scale as 


L(R e )=L(R 0 )[R e /R 0 V K2 3f) 

( 1 ) 

M(R e )=M (R 0 )[R e /R 0 ] fl(2 ~ 3f) 


Here Rq denotes Neptune’s mean orbital radius. The problem, therefore, is to find a 
mechanism that may lead to a low value for / and ensure uniform rotation. In addition, the 
mechanism should explain why the material is shed discretely, as individual gas rings, 
rather than as a continuous nebula spread throughout the equatorial plane. 
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In 1973, Prentice pointed out that if the interior of the proto-solar cloud (hereafter 
PSC) was very strongly convective such that in the cool outer layers the convective 
motions had mean radial speeds (v r ) that were many times larger than the local sound 
speed v_j = ^fyWfJp , then these motions create a turbulent radial stress in the gas, denoted 
by (^p t Vf 'j , which cannot be ignored. This stress can be much larger than the local gas 
pressure p gas = pv 2 jy = ptPT/ p . Here p, T, p are the local average gas density, 
temperature and molecular weight, p t is the mass density of moving material, 
y is the ratio of specific heats and £ is the universal gas constant. As the convective 
motions are driven by buoyancy, a simple scaling argument leads to the formula 

(p t vf) = /3pGM(r)/r (2) 

This expression is valid for a non-rotating cloud. Here M{r) denotes the mass interior to 
radius r, G is the gravitation constant, and p is a dimensionless constant called the 
turbulence parameter. The total radial stress at any point in the cloud is then given by 

Pm = p£T/p+PpGM{r)/r (3) 

At the center of the cloud, 0) = 0, T is a maximum and the turbulent stress is zero. At 

the surface M(r) achieves its maximum value, T is least and the ratio F t = {p t v t) /Pgas of 
turbulent stress to gas pressure is greatest. Thus F t increases steadily with increasing 
radius r. In other words, the inclusion of turbulent stress leaves the central region of the 
cloud unchanged but pushes the tenuous outer layers outwards. This causes the moment- 
of-inertia factor / to fall. Typically for a cloud of molecular H 2 , / drops from ~0. 1 1 2 to a 
plateau level of 0.010 as P increases from 0 to 0.1 (Prentice, 1973). When rotation of the 
cloud is included,/ falls by a further factor of (l + © e / 2) 2 , where ® e = co e Re /gM c 
denotes the ratio of centrifugal force to gravitational force at the equator and co e is the 
angular frequency. 

Thus, the inclusion of turbulent stress dramatically lowers the cloud’s moment-of- 
inertia factor, so accounting for the low masses of the planetary/regular-satellite systems 
relative to their primary bodies, as quantified in equation (1) and also equation (4) below. 
Prentice (1978a, 1978b) has also demonstrated that (i) thermally driven convective motions 
create a turbulent viscosity rft in the gas. This viscosity is essential for the cloud to rotate 
nearly uniformly and so efficiently convey angular momentum from the centre to its edge. 
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And (ii), that the breakdown of convective transport above the photosurface of the cloud 
leads to a scaling down of the turbulent stress there and the onset of a steep density 
inversion where the total pressure remains nearly constant. The highly turbulent interior of 
the cloud is thus capped with a dense shell of non-turbulent gas whose density exceeds that 
of the gas at the photosurface level by a factor (1+F,) - 30. 


The net result of the density inversion above the photosurface is that the turbulent 
cloud is able to store far more material at its edge than the non-turbulent structure. As the 
proto-solar cloud contracts inwards through the orbital radius R e = Rq of Neptune, where 
e e (* 0 )=i. a dense equatorial belt of non-turbulent gas is left behind, so forming the first 
orbiting gas ring. Its angular velocity C0q at radius Rq matches the local Keplerian 
frequency <a Kep =^GM 0 / Rq ■ The subsequent process of pressure restoration in the outer 
layers of the cloud then causes 0 e to drop catastrophically to a new solution value 0j < 1 . 
The equator of the cloud re-establishes itself at a new finite radius /?{ =y/?o(l + ©i/2). 

The contracting proto-solar cloud thus abandons a discrete family of orbiting gas rings at 
orbital radii R n (n = 0, 1, 2,. . .) which satisfy the equations 


K/R 


n + 1 


1 + - m " 


M n fn 


(4) 


Here m n denotes the ring mass and M n the residual cloud mass after detachment of the nth 
gas ring. If the contraction of the cloud mass occurs homologously, so that m n /M n and /„ 
stay reasonably constant, then the R n form a geometric sequence akin to the Titius-Bode 
law of planetary distances. 


We now apply equation (4) to the system of Galilean satellites, for which the 
observed mean geometric spacing ratio is (R n /R n +i) = 1-647. Setting M n = Mj, 
where M j = (l .89861 ± 0.00024)x 10 27 kg is Jupiter’s mass (Anderson et al., 1996a) and 
taking (f n ) = 0.02 , we obtain a mean gas ring mass (m n ) =1.075 xlO 25 kg . Assume, for 
the moment, that the gas consists of H 2 , He, and other elements in solar proportions, as 
compiled by Anders and Grevesse (1989) and Grevesse et al. (1992) and discussed in more 
detail below in Section 3. In that case, the total mass of condensable rock and H 2 0 ice in 
each ring is 
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"hock ~ 5.3 XlO 22 kg 
m H 2 0 =8.3xl0 22 kg 

Hence if the system of Galilean satellites did condense from rings of gas of roughly equal 
mass — dictated by the need for the proto-Jovian cloud to shed excess spin angular 
momentum during its inward radial contraction from the orbit of Callisto — we should 
expect to see rocky satellites of mass ~5xl0 22 kg and rock + ice satellites of mass 
-13x10 kg . Quite remarkably, the mass of the mostly rocky Europa is about this 
amount, namely 4.8 x 10 22 kg , and the mean of the masses of the two, rock+ice moons, 
Ganymede and Callisto, is 12.8 x 10 22 kg . 

It is extraordinary, therefore, how successful the modem Laplacian theory of solar 
system origin is able to account quantitatively even at the simplest level for both the broad 
orbital structure and distribution of mass in the Jovian system. This very satisfying feature 
was first noted by Prentice and ter Haar (1979). 


3. Gravitational Contraction and Thermochemical Development 
of the Proto-Solar Cloud 

3.1 Rate of radial contraction 

Consider now the physical structure and thermochemical evolution of the PSC during 
its gravitational contraction from an estimated initial radius /?,- = 7500/?©. This value 
exceeds Neptune’s mean orbital radius R 0 = 6463.2/?© (/?© = 6.9598 x 10 5 km). It is 
assumed that for equatorial radii R e < /?, , the cloud is sufficiently opaque that a convective 
equilibrium is established throughout its interior. The formation of a small central core of 
stellar density and mass -0.044M© is needed to release sufficient energy so that the levels 
of dissociation of H 2 and ionization of H and He throughout the rest of the cloud match 
the local thermodynamical equilibrium values (Prentice, 1978b). Here A/© = 1.98892 x 
10 30 kg is the solar mass (Standish, 1990), taking G = 6.67259 x 10' 20 km 3 s' 2 kg' 12 
(Cohen and Taylor, 1987). If the existence of the central core is ignored, then a full 
Kelvin-Helmholtz equilibrium cannot be achieved until the cloud radius R e has shrunk to 
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the value R* = 202.2 Rq, for the values of the other parameters given below. The rate of 
radial contraction for radius R e < R* is determined from the equation 
J = 2.8 1 13cr s RgTg = - d£ tot / dt . Here o s is the Stefan-Boltzmann constant and J 
denotes the luminosity of the fully rotating PSC whose shape is lenticular, with polar 
radius R p =^R e . Also, <f tot denotes the cloud’s total energy, including both the 
gravitational potential energy and dissociation/ionization energies. At radius R *, 
dR e /dt = - 2.1 66 x 1 0“ 5 v Kep , where v Kep = JgmTJI^ is the local Keplerian 
velocity for cloud mass M* = 1.0935M©. Since T e <*= R^?^ and, fairly closely, 
d <f tot / dt « (gM g / Rg )dR e I dt , it follows that 

x e (R e ) * 0.00371(R e /R Q ) OA /?©/ yr 

The cloud thus takes 8.3x 10 4 yr to contract from size /?, to R*. The contraction 
time from the orbit of Neptune to Mercury is 8.0 xlO 4 yr. In contrast, the contraction from 
Mercury’s orbit to radius R e = 1.5 Rq is much slower, taking ~ 10 7 yr . 

3.2 Angular momentum evolution and mass loss 

Next, we consider the physical structure and angular momentum evolution of the 
PSC . These are controlled by four parameters. The first of these is the initial cloud mass 
M t - prior to the shedding of the first gas ring, supposed at the orbit of Neptune. This is 
chosen so that the final cloud mass, at orbital radius R e = 1.5 Rq, is M©. The second 
parameter is the polytropic index ti\ of the super-adiabatic convective core of the PSC in 
which H is wholly atomic. Suppose n\ is treated as a variable and allowed to fall below the 
adiabatic value 1.5 and that the polytropic index «2 of the outer molecular layer is reduced 
by the same amount also. Then the gas pressure p n on the mean orbit of the gas ring shed 
at the orbit of Mercury increases relative to that for a purely adiabatic cloud. For gas ring 
temperatures in the vicinity of 1640 K, this increase in pressure increases the separation 
between the condensation temperatures of Fe-Ni alloy and MgSi 03 -Mg 2 Si 04 silicates 
(Grossman, 1972; Lewis, 1972). As a consequence, the condensate is starved of silicates 
relative to metals and so has a higher mean density (Prentice, 1991, 1995). Setting 
n t = 1.1, the mean density of the condensate matches the observed uncompressed mean 
density of Mercury, namely -5.3 g/cm (Anderson et al., 1987). 

The next parameter is the turbulence parameter (5. A value /?> 0.05 ensures that 
the cloud sheds discrete gas rings during its gravitational contraction. The choice 
f 3 = 0.08475 causes the peak in the condensate mean density due to the metal/silicate 
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fractionation process to occur exactly at Mercury’s orbit. Lastly, there is the ratio 
0 ph = [r ph / T c ]/ [u ph / // c ] of the photo-surface temperature to central temperature, divided 

by the corresponding ratio of molecular weights. Setting 6^ = 0.0027703 , the mean 

ratio (R n /R n+ \) of the orbital radii of adjacent gas rings between the orbits of Neptune to 
Mercury matches the observed mean value, namely 1.7231. For the above choices of n\ t (3 
and 0ph, we obtain M, = 1 . 1 69 M 0 . The PSC thus loses some ~1 5% of its mass during its 
contraction from Neptune’s orbit to present solar size. 

3.3 Methane and carbon synthesis in the outer layers 

During the homologous contraction of the PSC , the temperature at the equator scales 
closely as T e <x R~ 0 ' 9 . The temperatures T n of the gas rings shed at the orbits of Neptune, 
Uranus, Saturn and Jupiter are shown in Table 1 . The corresponding gas p n on the mean 
orbit of each ring are also shown. Within the outer layers of the PSC, the temperature and 
gas pressure increase with depth according to the super-adiabatic law p oc 7 2,9426 . Thus 
although at the surface of the cloud it is far too cold for any synthesis of methane CH 4 and 
C(s) solid carbon from the initial store of CO, this is not the case within the warm interior. 


Table 1: Physical Characteristics and Carbon Distributions of the Outer 

Proto-Solar Gas Rings 


Orbit 

Radial 

distance 

R„/Rq 

Elapsed 

time 

/io 4 yr 

Gas ring 
temperature 

T n /K 

Gas ring 
pressure 

p„ /bar 

Neptune 

6463 

0.80 

33.7 

8.48 x IO' 10 

Uranus 

4125 

2.81 

49.9 

5.72 x 10-9 

Saturn 

2050 

5.05 

92.4 

1.16 x IO' 7 

Jupiter 

1118 

6.36 

158.5 

1.64 x 10-6 


Carbon number fractions 


/ch 4 

/C(s) 

fc o 2 

f( CO 

Neptune 

0.1182 

0.0435 

0.2256 

0.6127 

Uranus 

0.3905 

0.1081 

0.1825 

0.3189 

Saturn 

0.7184 

0.1089 

0.0771 

0.0956 

Jupiter 

0.9093 

0.0367 

0.0255 

0.0285 
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As a consequence and as explained in detail elsewhere (Prentice, 1993), at sufficient depth 
inside the cloud CH 4 and, in some instances C(s), can be produced from CO and C0 2 
through catalytic synthesis on the surfaces of Fe-Ni grains. Some 49.5% of the Fe is 
available as catalyst at the assumed nominal synthesis temperature T cat = 475 K, after 
allowance is made for FeS formation. The reaction products are then assumed to be 
advected to the surface unchanged. It is assumed that the rapid inter-conversion between 
CO and C0 2 continues down to a quench temperature T q = T cat -50 K. The production 
rate of NH 3 from N 2 is negligible at temperature T cat for the low gas pressures of the outer 
layers of the PSC. (c/. Lewis and Prinn, 1980). As a consequence, all N remains as N 2 . 

Table 1 lists the distribution of carbon number in the four gas rings shed by the PSC 
at the orbits of Neptune, Uranus, Saturn and Jupiter. The contraction time T n from the 
initial cloud radius /?, = 7500 Rq to orbital radius R n is also given. The main point to note 
from this table is the steady conversion of CO and C0 2 to CH 4 in the outer layers of the 
PSC during cloud contraction. By the time the cloud has shrunk to the orbit of Jupiter, 

91% of the carbon has been converted to CFl 4 . Another point to note is the production of 
C(s). This is a transient phenomenon and only occurs for large cloud radii, where the outer 
layer gas pressures are very low. By the time the cloud has shrunk to equatorial radius 
R e ~ 870 Rq, all C(s) has gone. The condensate which forms at the orbit of the main-belt 
asteroids, i.e., at radius 605 Rq, is thus expected to be free of carbon (and also water ice). 
That is, the original population of planetesimals which condensed at the mean orbit of the 
asteroids was entirely stony. Inward orbital diffusion of carbon-bearing planetesimals at 
Jupiter’s orbit is needed account for the carbonaceous character of the outer asteroid 
population. 


4. Formation of the Pro to- Jovian Cloud 

4.1 Condensation of the proto-solar gas ring at Jupiter’s orbit 

Consider now the gas ring of mass 1.56xl0 28 kg that is shed by the PSC at Jupiter s 
orbit. It consists chiefly of H 2 and He. All other elements are present in solar proportion 
(Anders and Grevesse, 1989; Grevesse et al., 1992). Knowing the elemental abundances, 
we can construct the standard solar chemical mix, in which the elements are sorted into 
broad composition categories. This is given in Table 2. The total rock mass fraction is 
-^rock,© = 0.004968. 
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Table 2: Standard Solar Chemical Mix* 


Broad Category 

Species 

Mass fraction 

Gases: 

h 2 

0.7143498 


He 

0.2660000 


Ne, Ar 

0.0018185 

Ices: 

h 2 o 

7.6877 x 10' 3 


nh 3 

1.1297 x 10' 3 


ch 4 

4.0459 x 10' 3 

Rock oxides : 

MgO, A1 2 0 3 , Si0 2 , Fe 3 0 4 , 
CaO, Ti0 2 , Na 2 0, K 2 0, 
Cr 2 0 3 , P 2 0 5 , V 2 0 5 

3.7984 x 10’ 3 

Rock sulphides: 

FeS, NiS, MnS, ZnS 

1.1560 x 10' 3 

Rock halides: 

NaCl, CaF 2 

8.63 x 10' 6 

Unreacted metals: 

Co, Cu 

4.22 x 10' 6 

Rest: 


1.1 x 10' 6 


Based on: Anders and Grevesse (1989), Grevesse etal. (1992) and the 1987 Atomic Weights. 


It is now assumed that chemical condensation into solid grain occurs in accordance 
with the equilibrium condensation models first proposed by Larimer (1967), Lewis (1972) 
and Grossman (1972). A computational code has been specifically developed to deal with 
the gas ring configuration of the modem Laplacian theory of solar system origin (Prentice, 
1978a, 1990, 1991, 1995). At Jupiter’s orbit, the temperature of the gas ring at the moment 
of detachment from the PSC is 158.5 and the gas pressure on the mean orbit is 1.64xl0 - ^ 
bar. For these parameters, the condensate consists'simply of rock (mass fraction 0.6519), 
graphite (0.0146) and water ice (0.3335). The rock consists of MgSi0 3 (rock mass fraction 
0.2568), Mg 2 Si0 4 (0.1800), Si0 2 (0.0579), Ca 2 MgSi 2 0 7 (0.0423), MgAl 2 0 4 (0.0309), 

FeS (0.2044), Fe30 4 (0. 1755), and NiS (0.0229). The balance of the rock, making up mass 
fraction 0.0293, consists of NaOH, KOH, G2O3, P 2 C>5, Ti0 2 , V 2 C>5, MnS, ZnS, NaCl, 
CaF 2 , Co, Cu and residuals). 


Comparing the rock and ice mass fractions above with those of the standard solar 
mix, it is clear that the condensate at Jupiter’s orbit is water deficient by a factor of about 3 



Comparing the rock and ice mass fractions above with those of the standard solar 
mix, it is clear that the condensate at Jupiter’s orbit is water deficient by a factor of about 3 
relative to solar abundance. The reason for this shortfall is that the temperature of this 
proto-solar gas ring is just 5.1 K less than the condensation temperature of H 2 0 on the 
mean orbit. Thus only 33.5% of the available water vapour of the gas ring can condense as 
ice. Allowing for the formation of NaOH and KOH, the net condensed mass of H 2 0 is a 
fraction ^ H; , 0 = 0.3495 of the initial free H 2 0 vapour mass. Hence the very same process 
of chemical fractionation which operated in the gas ring at Mercury’s orbit, causing a 
shortfall in the condensed mass fraction of silicates relative to metal there, also operates 
with the same effect on the ice and rock mass fractions at the orbit of Jupiter. 

4.2 Formation of a planetary core 

After chemical condensation into solid grains has occurred, the grains migrate from 
all regions of the residual gas ring onto its mean orbit to form a concentrated torus. This 
process is depicted in Figure 2. Self-gravitation within the grain torus results in the 
formation of planetisimals which remain bound to the mean orbit through the combined 
action of gas drag and the gravity of the PSC. These processes have been examined 
elsewhere (Prentice 1978a, 1980; Hourigan, 1997). The process of accumulation of the 
first generation of planetisimals to a single orbiting planetary core remains a largely 
unexplored aspect of the modem Laplacian theory. 

The total mass of condensate at Jupiter’s orbit is 19.9 M®. Assume now, in line with 
the calculations of Mizuno (1980) and others, that the first stage of Jupiter’s formation was 
the accumulation of a solid planetary core of mass roughly 15 M®. This core then serves 
as sink for capturing the light gases and residual solids of the gas ring at Jupiter s orbit. 
These solids are easily captured since they form a concentrated slurry along the mean orbit. 
In contrast, probably only a small fraction of the gas can be cap-tured because of radial 
dispersal due to viscous forces induced by the steep gradient of the orbital angular velocity 
distribution in the initial gas ring. We have ffl(5,z)=©o(^o/ i ) 2 • Here (s,z) refer to 
cylindrical polar coordinates, s = Rq, and z = 0 defines the mean orbit, as shown in 
Figure 2, and a> 0 is the angular velocity on the mean orbit. 
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Figure 2. A simple schematic view of a gas ring orbiting the protosun. As the gas cools, solid grains 
condense out and migrate to the mean orbit of the ring under the combined influence of the 
protosolar gravity and drag induced by motion through the gas. The mass of gas is much 
less than the mass M of the protosun and it has uniform specific orbital angular momentum. 

The orbital angular velocity co thus varies as coocl/s 2 , where s is the cylindrical distance 

from the common axis of rotation of the system. On the mean orbit S = R 0 , co equals the 

local Keplerian value <o 0 = a Ktp = Jgm / r\ ■ 

4.3 Composition of the primitive Jovian atmosphere 

Hence the primitive Jupiter probably consisted of a heavy-element core of mass 
m core ~ 1 53/© surrounded by a radially extended envelope of mass m env » 3203/©, 
containing ~53/© of the condensate mix defined above, along with the uncondensed 
chemical constituents of the proto-solar gas ring, namely H 2 , He, Ne, Ar, N 2 and CH 4 , and 
including the remaining H 2 0 vapour of mass ~1 .63/®. This data imply a rock mass 
fraction Z rock « 0.010 and total H 2 0 mass fraction Z H2 q a 0.010. Comparing these 
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fractions with the solar abundances given in Table 2, we see that the rock is enhanced by a 
factor of 2 and H 2 0 by a factor of 1 .3, relative to solar expectations. 

On the basis of the above model for Jupiter’s formation, Prentice (1996a) predicted 
that the sulphur content of the Jovian atmosphere to be twice solar, yielding a H 2 S to H 2 
number ratio g / ^h 2 = 7.4x 10 . This prediction was subsequently confirmed by the 

Galileo Probe. In their initial report, Niemann et al. (1996) found that [«h 2 s/ 77 h 2 1 J = 

(7.7 ± 0.5)x 10 -5 . This result is a very significant triumph for the modem Laplacian 
theory since H 2 S is a true indicator of the heavy element abundance of the Jovian 
atmosphere. The H 2 S measurement has been re-confirmed by subsequent re-analysis of 
the Probe data (Niemann et al., 1998). Although the Galileo Probe measurements were 
unable to give any estimate of the global H 2 0 abundance, the predicted enhancement of 
1.3 times solar abundance made by Prentice (1996a) is consistent with the 1-2 times solar 
value deduced indirectly from meteorological considerations (Showman and Ingersoll, 
1998). The predicted near solar abundance for He has also been confirmed by the Probe 
data (Zahn et al, 1998), but the -3 times solar abundance of CH 4 still remains 
unexplained. 

Lastly, in quantifying the chemical composition of the proto-Jovian envelope it is 
convenient to let C en ~ 2 denote the enhancement in abundance of fully condensed species 
relative to solar abundances in the material captured from proto-solar orbit. Rock is thus 
enriched by factor Z en , while the unnormalized H 2 0 abundance is given by 

2fj(H 2 0) = (l + (^ en -1 ¥h 2 0 k© (H 2 0) (5) 

where ^ HO * 03495 is the fraction of condensed H 2 0 in the proto-solar gas ring. 
Choosing £ en = 2.0812 , the H 2 S to H 2 number ratio in the proto-Jovian envelope exactly 
matches the Probe measurement. Table 3 gives the predicted broad chemical composition 
of the proto-Jovian atmosphere for this choice of C en . Because of the high gas pressures 
in this atmosphere compared to those of the PSC, it is safe to assume that all incoming 
C(s), CO and C0 2 are converted to CH 4 , and all N 2 to NH 3 (Lewis and Prinn, 1980). The 
abundances of both C and N in the proto-Jovian atmosphere are again expected to be 
nearly solar. The H 2 0 abundance of the Jovian atmosphere is now 1.35 times solar for this 
choice of <^ en . 
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Table 3: Proposed Broad Chemical Composition of the Proto-Jovian 
Atmosphere: £ en = 2.0812 


Broad category 

Species 

Mass fraction 

Gases 

h 2 

0.70846 


He 

0.26382 


Ne, At 

0.00180 

Ices 

h 2 o 

1.0373 x 1 O' 2 


nh 3 

1.121 x 10' 3 


ch 4 

4.172 x 10‘ 3 

Rocks 

Same as Table 2 

1.0255 x 10* 2 


5. Gravitational Contraction of the Proto-Jovian Cloud 

We are now in a position to model the contraction of the proto-Jovian cloud 
(hereafter PJC). Since the orbit of the outermost regular satellite Callisto is 26.34 R } , 
where /?j = 71 ,492 km is the 1-bar equatorial radius of Jupiter (Davies et al., 1992), it is 
reasonable to assume an initial equatorial radius for the PJC, R e = ~ 50Rj. The cloud is 

assumed to acquire a super-adiabatic convective equilibrium at that radius through heat 
generated by gravitational contraction. We also assume that the orbit of Callisto marks the 
radius R 0 where the centrifugal force at the equator of the PJC first matches the 
gravitational force and a gas ring is shed. Thereafter, during subsequent contraction, we 
assume that gas rings are shed close to the present mean orbital radii of the other large 
satellites at the time of detachment from the cloud. Next, in computing the structure of the 
PJC, we adopt the same values for the parameters (5 and n.\ that were derived for the PSC, 
namely /? = 0.08475 and n i = 1.1. The parameter 0^ = 0.0049066 is chosen, as before, so 
that the mean geometric ring spacing ratio (R n /R n+ \) matches the observed value, viz. 

1 .6465. The initial cloud mass is M i = 1 ,0542Mj. The contraction occurs nearly 
homologously from radius to 3% A steady but modest departure from uniformity 
occurs through the progressive dissociation of H 2 in the cloud interior. The temperature 
7p h at the photosurface increases closely as T ph « R^ 9 , ignoring solar radiation. 
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Allowing for heating due to the primitive sun, whose estimated luminosity is ~ 0.7 =/©, 
the background temperature at Jupiter’s orbit for an assumed albedo of 0.5 is 7b = 93.8 K. 
The effective mean temperature of a gas ring at the moment of detachment from the PJC is 
then 

t„ = T ctt = ().5r p 4 h + r b 4 ) 0 - 25 (6) 

The factor of 0.5 multiplying 7^ in equation (6) is a geometric factor associated with the 
fact that the surface area for radiance from a gas ring of given minor radius is twice the 
area receiving energy from the PJC. 

For cloud radii R e = 3i?j, T t ff = 898 K. Since further homologous contraction 
results in unreasonably large T e ff it is assumed that for R e < 3Rj, begins to decline in a 
manner such that the cloud’s moment-of-inertia factor stays constant, namely/= 0.03875 
while T e ff smoothly passes to a final cloud temperature r end = 750 K, corresponding to 
radius R e = 1.5/?j. The solid heavy curve in Figure 3 gives the plot of the effective 
blackbody temperature Tgff at the equator of the PJC as a function of the cloud equatorial 
radius R eq = R e . The open circles mark the temperatures of the gas rings cast off at the 
orbits of the four Galilean satellites. The lighter curves define the condensation 
temperatures of the major chemical constituents on the mean orbit of the locally detached 
gas rings, plotted against cloud radius R eq . These curves correspond to the choice £ en = 2. 
They are barely distinguishable from those for the case ^ en = 2.0812 which is adopted for 
the rest of the paper. 

Table 4 gives the gas ring temperatures T n , mean orbit pressures p n and detachment 
times for the four Galilean gas rings, along with the broad chemical composition of the 
condensate and both the total mass of condensate and its mean specific density p s . This 
last quantity is calculated for temperature 298.15 K and 1 bar pressure. 
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Table 4: Broad Characteristics of the Family of Proto-Jovian Gas Rings 


Orbit 

Radial Elapsed Gas ring Gas ring 

distance time temperature pressure 

R n /R J f„/10 6 yr T n / K p n i bar 

Io 

Europa 

Ganymede 

Callisto 

5.900 6.46 463 18.3 

9.387 5.31 301 2.77 

14.973 3.76 199 0.405 

26.335 1.44 129 0.0399 


Condensate mass fraction 

Condensate Specific 

mass density 

(102 6 g) p s /gcnr 3 

Rock H 2 0 NH 3 

Io 

Europa 

Ganymede 

Callisto 

1.000 0 0 

1.000 0 0 

0.5249 0.4751 0 

0.4956 0.4569 0.0475 

1.76 3.719 

1.88 3.134 

3.37 1.477 

3.33 1.418 


6. Chemical Condensation of the Galilean Satellites 

6.1 The general picture 

Inspecting Figure 3 and Table 4 we observe that the gas ring model of planet/satellite 
formation does a very good job of explaining the broad chemical structure of the Galilean 
satellite system. As one passes from the orbit of Io to Callisto, the mean density of the 
condensate declines with increasing orbital radius in a manner that closely replicates the 
observed trend in satellite mean density, which is shown in Table 5. The agreement is even 
better after allowance is made for the effects of self-compression (see section 8). The 
condensation point of water ice lies about midway between the orbits of Europa and 
Ganymede. This explains the sharp drop in mean density between the inner and outer pairs 
of satellites. The condensation model thus suggests that Europa possessed no initial store 
of H 2 0 ice, while Ganymede and Callisto acquired roughly equal proportions of rock and 
ice. In addition, Callisto contains some 5% by mass of NH3 incorporated in the 
condensate as the mono-hydrate ice NH3 H2O. 
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Figure 3. Effective blackbody temperature T c ff at the equator of the proto-Jovian cloud (PJC) versus 
cloud radius R e q which is expressed in units of Jupiter's present equatorial radius 
Rj =71,492 km . The temperatures of the gas rings at the present orbits of the Galilean 
satellites are shown with the open circles. The other lines in the figure give the condensation 
temperatures of the principal rocks and ices. These are computed for the gas pressure on the 
mean orbit of the locally detached gas ring 


It is remarkable that the same choice of turbulence parameter/? which accounts for 
the broad chemical structure of the planetary system, also accounts for that of the Galilean 
system. There is, however, one problem with Table 4 in that the predicted mass of 
condensate in each gas ring is about twice the observed satellite mass. These are listed in 
Table 5. This means that either the efficiency of the gravitational accumulation process 
was only 50% or, more likely, that the real structure of the proto-Jovian cloud and its 
family of gas rings is not completely represented by the simple model for the radial 
turbulent stress, expressed in equation (2). That is, the gas ring masses may be 
overestimated by a factor of two. An attempt to improve the theory of supersonic turbulent 
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convection especially to model the degradation of turbulent stress in the outer layers of the 
cloud is being currently undertaken. This study takes into account the results of a detailed 
numerical simulation of supersonic turbulent convection in a model solar atmosphere (Dyt, 
1997; Dyt and Prentice, 1998). 


Table 5: Bulk Physical properties of the Galilean Satellites 


Satellite 

Mean physical radius 3 (km) 

Mass b ’ c (10 22 kg) 

Mean density(g/cm 3 ) 

Io 

1821.3 ±0.2 

8.9319 ±0.0012 

3.529 ± 0.001 

Europa 

1565 ± 8 

4.8000 ± 0.0006 

2.990 ± 0.046 

Ganymede 

2634 ± 10 

14.8167 ± 0.0020 

1.936 ±0.022 

Callisto 

2403 ± 5 

10.762 ±0.005 

1.852 ±0.012 


a Radii from Davies et al. (1992) 

b Masses: Io (Anderson et al. 1996b), Europa (Anderson et al., 1997a), Ganymede (Anderson et al., 
1996a), Callisto (Campbell and Synnott, 1985) 

c Gravitational constant: G = (6.67259 ± 0.00085) x 10' 20 km 3 s' 2 kg' 1 (Cohen and Taylor, 1987) 


6.2 Iron-rich Io 

Table 6 gives a breakdown of the predicted rock constituents of each of the Galilean 
satellites. Io rock differs markedly from that of the other three satellites, which are all very 
similar in composition. It is much drier, containing only some 0.5% by weight H9O tied up 
in tremolite. A second important difference with Io lies in the distribution of iron. All of 
the satellites have a common supply of FeS which accounts for roughly 50% of the Fe 
atoms. For the outer three moons, the remaining iron is all present as Fe 3 0 4 . For Io, the 
gas ring is too hot for Fe to be rusted by water vapor to form Fe 3 0 4 . Instead, fayalite 
(Fe 2 Si0 4 ) appears forms via the solid-state diffusion reaction. 

2Fe(s) + Si0 2 (s) + 2H 2 0(g) = Fe 2 Si0 4 (s) + 2H 2 (g) 

This has a formation temperature 7f orm = 510 K in the Io gas ring. An alternative reaction 
path, involving the diffusion of Fe with MgSi0 3 , cannot proceed here since it has a 
formation temperature of 420 K ( cf. Fegley, 1988). Hence there is insufficient Si0 2 (s) left 
after the much swifter tremolite reaction (T trem = 1000 K) to use up all the remaining Fe. 
Allowing for FeTi0 3 formation as well, some 16.5% of the initial Fe remains as free metal 
in Io. 
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Table 6: Predicted Initial Rock Constituents of the Galilean Satellites 



Mass fractions X t 

Chemical Species 

Io 

Europa 

Ganymede 

Callisto 

Ca 2 MgSi 2 07 

0 . 

0.0389 

0.0387 

0.0387 

MgAl 2 0 4 

0.0317 

0.0284 

0.0283 

0.0283 

MgSi0 3 

0.1997 

0.0003 

0 . 

0 . 

Mg 2 Si0 4 

0.1848 

0.0002 

0 . 

0 . 

Si0 2 

0 . 

0.2622 

0.2615 

0.2615 

Ca 2 Mg 5 [Si 8 0 22 ] (OH) 2 

0.1293 

0 . 

0 . 

0 . 

Mg(OH) 2 

0 . 

0.2739 

0.2733 

0.2733 

NaOH + KOH 

0.0121 

0.0108 

0.0108 

0.0108 

Fe 

0.0434 

0 . 

0 . 

0 . 

FeS 

0.2366 

0.2068 

0.1911 

0.1911 

Fe 2 Si0 4 

0.1264 

0 . 

0 . 

0 . 

Fe 3 0 4 

0 . 

0.1448 

0.1574 

0.1574 

Ni 

0.0152 

0.0100 

0 . 

0 . 

NiS 

0 . 

0.0056 

0.0209 

0.0209 

Residuals 

0.0208 

0.0181 

0.0180 

0.0180 

TOTAL 

1.0000 

1.0000 

1.0000 

1.0000 

Mean density (g/cm 3 ) 

3.7192 

3.1342 

3.1334 

3.1334 


6.3 Europa, Ganymede and Callisto 

Consider now Europa, Ganymede and Callisto for which it is argued, as follows, that 
initially all Mg is tied up chiefly as brucite Mg(OH) 2 , rather than MgSi03-Mg 2 Si0 4 . 
Within a rising convective updraft inside the proto-Jovian cloud, both MgO(s) and 
MgSi0 3 (s) begin to condense as solid grains as soon as the temperature falls below 
-1950 K. MgO(s) is likely to be the dominant condensation product, however, because of 
the fewer molecular collisions required. Next, Si0 2 (s) does not begin to condense until the 
updraft has cooled to 1 850 K. But the time taken for an updraft travelling with the local 
sound speed = yfyp ] p to pass from the 1950 K to 1850 K radius levels, when the PJC 
has Europa’ s orbital dimension, is only 1250 s. Thus whilst it is reasonable to assume that 
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all vapor-vapour reactions, such as Mg(g) + H 2 0(g) — » MgO(s) + H 2 (g), pass to 
completion in that interval, solid-vapour reactions, such as the reaction 
MgO(s) + SiO(s) + H 2 0(g) — > MgSi0 3 (s) + H 2 (g) may not have time to proceed far 
before the SiO vapor condenses to Si0 2 (s). This is because the time-scale for solid- vapour 
reactions exceed that for vapor- vapour reactions by the collision factor A^. 3 - 10 3 -10 4 , 
where is the number of atoms per grain (Fegley, 1988). depends on the grain size 
a gr Typically ~ (0. l-l)//m . Hence a large portion of the Mg advects to the cloud 
surface as MgO(s). Since the formation temperature for brucite at Europa’s orbit is 414 K 
(see Figure 3), the conditions are favorable in the outermost cloud layers for the copious 
hydration of MgO grains. This solid-vapour reaction should proceed far more swiftly them 
the solid-state reaction MgO(s) + Si0 2 (s) — > MgSi0 3 (s) which requires intimate grain 
contact for long times whilst solid state diffusion takes place. 

Returning briefly to Table 6, the last row gives the mean density of the rock 
condensate calculated at 298.15 K and 1 bar pressure. Amongst the residual rock 
constituents in each satellite, in order of decreasing abundance, are MnSi0 3 , Cr 2 0 3 , P 2 0 5 
and NaCl. This last mineral makes up 0. 16% by mass of Io rock and 0. 14% of the other 
satellites’ rock. 


7. Construction of Thermal and Structural Models for the Satellites 

The final stage of this investigation is the construction of detailed thermal 
evolutionary and structural physical models for each of the Galilean satellites. To this end, 
two large computational codes have been developed over the past several years in readiness 
for the Galileo mission. 

7.1 The thermal evolutionary code 

The first of the computational codes determines the temporal evolution of the 
temperature distribution within both homogeneous and differentiated satellites of arbitrary 
chemical composition. The only heat source considered is radiogenic energy released 
through the decay of K 40 , Th 232 , U 235 and U 238 . Io is thus excluded from this study since 
its heating is derived mostly from tidal action with Jupiter. For simplicity, all satellites are 
assumed to be solid and, initially, to have a uniform temperature T s throughout, equal to 
the present day surface value T s = 103 K (Lupo and Lewis, 1979). The mass and radius of 
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each model are chosen to match the present day values (see Table 5). Empirical formula 
have been constructed to represent the temperature dependence of the specific heats 
Cp ■ (T) and thermal conductivities K t (7) of the major rock-bearing minerals listed in 
Table 6, and those of the ices H2O and NH3. Here i is the species index. With respect to 
the minerals, the data compilations of Touloukian et al. (1970 a,b) and Robie et al. (1978) 
are very useful. As a rule of thumb, thermal conductivity data is scarce, especially at low 
temperatures (T < 298.15 K). It appears non-existent for FeS, NiS, Ca 2 MgSi 2 0 7 , 
Mg(OH) 2 , NaOH, KOH, and the alkali feldspars NaAlSi 3 0 8 and KAlSi 3 0 8 . To get around 
this difficulty, a nominal conductivity is prescribed, given by 

K t (T) = 4.2 [pi (298.15)/p std ]k 2 (. T)/k 2 (298.15) Wm^K' 1 (7) 

-1 -1 'X 

This formula is based on the standard rock conductivities K std = 4.2 Wm" K' and specific 
density Pstd = 5 5 g/cm ~ adopted by Ellsworth and Schubert (1983). In addition, a 
temperature dependence on K[ is the imposed which is the same as that of MgAl 2 0 4 
(i.e., i = 2) (Slack, 1962; Kingery et al., 1954). 

Because H 2 0 ice is a major constituent of Europa, Ganymede and Callisto, it is essential to 
have accurate thermophysical data for this substance. Table 7 gives empirical 
representations for C p (7), k{T) as well as those of other properties of H 2 0 ice to be used 
later on. The data used for representing k(T) is in excellent agreement with the more 
recent low-temperature measurements of Ahmad and Phillips (1987). All of the data in 
Table 7 is for phase Ih of H 2 0 ice. Data for the other ice phases does not seem to exist. 


It is assumed that all heat is conducted outwards according to the heat diffusion 


equation 


pc 


■ d JA - 

p dt 


1 d 

r 2 dr 


r 2 K 


dT d 

dr 


\ 

) 


+ ^ nuclW 


( 8 ) 


Here T d ( r,t ) denotes the diffusive temperature at radius r and time t and p,C p {T d ) and 
jt"(7j) are the mean density, specific heat and thermal conductivity. For simplicity, the 
temperature dependence of p is ignored at this stage. It is calculated at the adopted 
reference temperature 7 ref =200K. Within the satellite, only the variation of p with 
changes in composition is taken into account. The volumetric radiogenic heat production 
rate H nuc \(t) is taken from Turcotte and Schubert (1982, p. 140). 
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Table 7: Thermophysical Properties of H 2 0 Ice (Phase Ih) 


Quantity 

Symbol 

[Unit] 

Temperature 
range [K] 

Formula 

Data source 

Cubical 

expansion 

coefficient 

a(D 

IK"'] 

(20, 100.15) 
(100.15, 273.15) 

-6.94 x 10-7 r+ 1.0818 x 10-8 7*2 

-3.700 x 10-5 + 7.79 x 10* 7 T- 2.022 * 10- 10 T 2 

Leadbettcr(1965) 

[See also Hobbs (1974), p. 350] 

Expansion 

coefficient 

integral 

/i(r)=£a(r>* 

(20, 100.15) 
(100.15, 273.15) 

-3.47 x 10- 7 r 2 + 3.606 x10-9 7*3 
7.8 x 10-6 - 3.700 x 10-5 ?+ 3.395 * jq-7 j2 
-6.74 x 10 -n T 2 


Adiabatic 

compressibility 

Aden 

[bar 1 ] 

(60, 110) 
(110, 133.15) 
(133.15, 273.15) 

9.1884 x 10-6 /(l -2.7157 x 10-4 T- 1.5933 x 1( ^6 
3.2196 x 10-6+ 5.015 x io-«r 
1.945 x 10-5/(1- 1.7305 x 10-3 r c - 1.804 x \<y6 
[f c = r- 273.15] 

Proctor (1966) 

[see also Hobbs (1974), p.263] 
Dantl (1968) 

[see also Hobbs (1974), p.261] 

Specific heat 

C p (7) 

[kJ/kg.K] 

(17, 85) 
(85, 170) 

(170, 273.15) 

-0.13779+ 1.32002 x 10*2 T- 3-23053 x ]0-5j2 
-0.78302 + 2.92186 x 10*2 T- 1.65222 x I(H T 2 
+ 3.97284 x 10-7 

0.12996 + 7.2419 x 10-3 j 

Giaque and Stout (1936) 

Thermal 

conductivity 

k(T) 

[W/m.K] 

(16, 100) 
(100, 273.15) 

1.605 + 305.3 /f + 9220 iT 2 
0.330 + 525 fT 

Klinger and Neumaier (1969) 
Dillard and Timmerhaus (1965) 


Next, T d ( r,t ) is compared at each time step with the local solid-state convective 
creep temperature 7" creep given by 

_ J0.6r me i t (p) [rock + ice mix] 
creep [1073.1 5K [rock only] (9) 

Here TVneU (?) denotes the melting temperature of the major ice component, namely HtO. 

It depends on the local pressure p = p{r) which is estimated from a simple integration of the 
hydrostatic support equation. If T d ( r,t ) rises above r creep at any time, we set T = T c 
The 0.6 factor in the creep criterion for ice is taken from Ellsworth and Schubert (1983). 
The rock creep temperature is from Schubert et al. (1979; cf Ransford et ai, 1981). 

Lastly, all satellites are evolved for 4600 Myr. The resulting temperature profile is used in 
the computation of the internal physical structure. 

7.2 The structural code 

The computational code for the satellite structure is devised as follows. First the 
model is divided into distinct compositional zones. Second, the pressure at the boundary 
between convective and wholly conductive ice zones in the thermal model is used to flag 
the corresponding boundary in the structural model. The ice zones are then further sub- 
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divided according to the pressure boundaries between the four relevant phases of water ice. 
These boundaries depend on the local temperature 7 and are given by the equations 


Plh/ll(T) — 84.7+8.6287 

bar 


PlWi(T) = 11325.-24.07 

bar 

(10) 

7vi/VIIlCn = 9636.+41.77 

bar 



These formulae have been derived by least square fits to the data of Bridgman (1912), 
Whalley (1969) and Brown and Whalley (1966), respectively. Knowing the temperature 
profile of 7 versus p within the thermal model of each satellite, equations (10) then yield a 
unique solution for each of /?ih/n, 7II/VI anc * Pvi/VlII- 

Once all the compositional and phase boundaries are specified, a uniform mean 
temperature T j is assigned to each zone j. In the case of conductive layers, Tj is esti- 
mated from the thermal code temperature profile. For convecting layers, Tj can be deter- 
mined by an iterative process within the structure code, by making use of equations (9). 

7.3 Computation of specific densities 

Lastly, the local density p =p (p,7) at each point in the satellite model is given by the 
equation 

yp(pj)=f,x t /Pi{ P j) 

i=i 

where X t and p ( (p,T) denote the mass fractions and specific density of species / and n 
denotes the total number of species. For rock species, the p, ( p,T) are estimated by the 
semi-empirical formulae 

pi{pj)= Pi{lT){\+ K' 0 (p-\)l K 0 ) U K '° (11) 

Here p ; (1,7) denotes the specific density of species i at temperature 7 and 1 bar pressure, 
and K 0 and K ' 0 are the bulk modulus and its pressure derivative at zero pressure and 
temperature 298.15 K . Space does not permit a full description of these quantities in this 
paper. However, several important points deserve mention here. First, with regard to 
Pl (1,7) excellent data is given for most rock minerals in Robie et al. (1978) for room 
temperature 7 0 = 298. 15 K. Second, Touloukian et al. (1970c) provide an excellent source 
of thermal expansion data to assist in the computation of p, (1,7) for 7 > 298. 15. With the 
exceptions of Fe, Ni, Fe 3 0 4 and Si0 2 , however, very little thermal expansion data exists 
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for T < 298. 1 5. For low temperatures, therefore, a simple extrapolation is assumed via the 
linear expansion formula 

( 12 ) 

Here (AZ) ; denotes the change from the room temperature length Lg at temperature T. The 
exponent 1.7 provides a best fit for the Fe thermal expansion data. The coefficient Bj is 
chosen so that (AL)JLq smoothly matches onto the high temperature data. 

Lastly, a useful compilation of Kg and Kg values for several minerals is given in 
Prentice (1996b). Data for Mg(OH) 2 is given in Fei and Mao (1993). Kusaba et al. (1997) 
report improved compressibility data for FeS, especially in regard to two phase transitions 
of this substance at pressures of 34 kbar and 67 kbar, respectively. Lupo and Lewis ( 1 979) 
were the first to provide detailed equations of state [EOS] for the different phases of H 2 0 
ice for use in satellite modelling. Much of the data they assembled, however, pre-dated the 
advent of X-ray diffraction measurements. Table 8 is an attempt to provide an improved 
EOS for water ice using the diffraction data. Because of the scarcity of thermal expansion 
data for the higher pressure phases of ice, it has been necessary to assume that the specific 
density, /?h 2 OJ (p^) is a separable function of pressure p and temperature T. Here j is 
the phase index. A temperature dependence is then imposed which is the same as that of 
ice Ih, as detailed in Table 7. 


Table 8: Equation of State of H z O Ice: Pa 2 oj (p,T) = Pn 2 o,j (M)e fl(T) Fj ( p ) 


Icc 

Phase 

Index j 

Density at I bar and 0 K 
j (1,0) g/cm 3 

Compression factor 
Fj{p) 

Data Source 

lh 

I 

0.93298* 

e As° (/’->)_ /? js0 (T ) + chT)T/ Pl (l,T)C p (T) 

Table 6 


2 

1.18066 

2 ~ 1-21736x10”^ 

Kamb (1964) 

Bridgman (1912) 
Landolt-Bomstein (1975, p.5) 

VI 

3 

1.31408 


1-7.5282 X 1 0 6 fr-l) + 4.068 x 1 <f ' ‘(p-l) 2 J 

Kamb (1965) 

Bridgman (1937) 
Landolt-Bomstein (1975, p.7) 


■ 

1.49084 


l-4.3194xl0~ 6 (p-l) 2 + 2.1036x1<f n fr-0 2 

Bridgman (1937) 

Hobbs (1974, p.79) 
Landolt-Bomstein (1975, p.7) 


♦This value chosen so that Ph 2 o.i( ] , 273.15) = g/cm 3 [Ginnings and Comiccini (1947); see also Hobbs (1974, p. 346)] 
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8. The Satellite Models 


8.1 Io 

The observed mean density of Io, namely 3.529 ± 0.001 g/cm- 5 (Anderson et al, (1996b) 
lies well below the uncompressed solid rock density expected at Io’s orbit according to 
Table 6. It would therefore appear that the bulk composition of Io cannot be accounted for 
by the condensation model presented in this paper. This result is not a surprise however 
since, as mentioned in Section 1, most scientists agree that at least part of the interior of 
this volcanically active moon is molten. Melting produces a slump in specific density of 
most rocks, particularly in FeS. Prentice (1996b) thus constructed a pre-Galileo model for 
Io that is consistent with the observed mean density. All FeS, including dissolved Ni, Cu 
and Co are accumulated at the center to form a molten core of mean density ~4.0 g/cm 3 
and mass -25% of the Io mass. This model yields a predicted mean axial moment-of- 
inertia coefficient for Io of/j 0 = 0.390 ± 0.002. 

Now the Galileo Orbiter observed f lo to be 0.378 ± (Anderson et al, 1996b). 

This is -1 Via shy of the predicted value. The real Io thus contains heavier material at its 
center — and correspondingly lighter material in its mantle — than that predicted in Prentice 
(1996b). It is hoped that the revised compositional suite for the initial Io presented in 
Table 6 may provide the basis of a 4-zone model for Io that is consistent with all of the 
Galileo gravity data. This model — yet to be constructed— consists of an innermost core of 
solid Fe-Ni-Co-V (of mass percent -6% of the Io mass M Io ) surrounded by a layer of 
molten FeS-Cu (mass -24%), then a layer of molten Fe 2 Si0 4 (mass -12%%) and lastly an 
outer, iron-free, mantle of calcium and magnesium silicates. 

8.2 Europa 

Figure 4 shows the internal temperature distribution T(r,f) in Europa (E) plotted as a 
function of radius r, expressed in units of the satellite radius at several chosen times t in 

the evolution. Starting from an initial uniformly cold state, the interior of the satellite 
rapidly warms under the influence of radioactive decay. After 500 Myr, the temperature in 
the core has already reached the value -800 K where decomposition of brucite Mg(OH) 2 
can take place (Ransford et al, 1981). We therefore assume that complete dehydration of 
the initial Europa rock takes place early in the evolution of the satellite, leading to a present 
day model consisting of a dehydrated rock core surrounded by an outer mantle of H 2 0, 
including dissolved salt. 
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Figure 4. Internal temperature distribution of a 2-zone model of Europa versus fractional radius x = r/R E 
for several key times during the course of its radiogenic evolution to present age (4600 Myr). 
The satellite model consists of a dehydrated rock core of radius Xj = 0.905 and uniform density 
3.68 g/cm 3 , capped by a solid H 2 O mantle. After time 50 Myr, the temperature at the rock/ice 
boundary rises to the value 0.6 T^ei, (pj) = 151 K where the ice becomes unstable towards 
solid-state convection. Here is the melting temperature of H-,0 ice and p t = 2.01 kbar is 
the boundary pressure. Thereafter the temperature at this boundary stays constant. The surface 
temperature of the model is set at 103 K throughout the evolution. 


If NaOH and KOH are converted to alkali feldspars and all NaCl is evicted with 
H 2 0, then the final mass of the dehydrated rock core is 0.91 16M E . The rock has the same 
relative composition as that listed in Table 6 except that the entries for MgSi0 3 , Mg 2 Si0 4 , 
MgAl 2 C>4 and Si0 2 change to 0.2334, 0.1766, 0.0097 and 0.0000, respectively. 

NaAlSi 3 Og and KAlSi 3 Og make up a fraction 0.0692 of the total satellite mass. The 
uncompressed mean density of the new rock at room temperature (298.15 K) and 1 bar 
pressure— hereafter RTP— is 3.6787 g/cm 3 . NaCl makes up 1.63% by mass of the outer 
H 2 0 ice mantle of total mass 0.0884 M E . 

The temperatures in the outer layers of the satellite generally peak around time 
-1500 Myr. The radius and pressure at the rock/ice boundary are rj = 1415 km and 
Pl =2110 bar. The temperature T(r,t) at that point very quickly rises to the value - 15 1 K 
where the ice becomes convective. Interestingly, this pressure is very close to the value 
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-2100 bar where the melting temperature 7 melt ip) of H 2 0 ice is a minimum. Thus as 
physical considerations dictate that 6T{r,t)/6r< 0 , it is necessary to replace the creep 
criterion in equation (9) by T aeep = 0.6 r melt [j> x ) in regions where 6 7 me it (p)/° P <0 - M 
time t = 4600 Myr, the thickness of the inner layer of convecting ice is - 76 km. This is 
just half of the total ice layer thickness. 

The structural parameters of the computed present-day Europa model are given in 
Table 9. Figure 5 is a simple schematic picture of the proposed satellite model. The mean 
density of the 2-zone rock/ice model presented here agrees well with the observed mean 
density p E = 2.990 ± 0.046 g/cm 3 , for the satellite radius R E = 1565 km derived from 
Voyager measurements (Davies et al., 1992). The moment-of-inertia coefficient is also in 
good agreement with the recommended value of 0.346 ± 0.005 obtained by the Galileo 
Orbiter gravity experiment (Anderson et al., 1998a). In fact, a value of 0.347 ± 0.004 
predicted prior to any of the Galileo flypasts (Prentice 1996c) coincides exactly with the 
measurement of 0.347 ± 0.014 obtained during the E6 flypast (Anderson et al., 1997a)! 
Thus a simple 2-zone Europa model, consisting of a dehydrated rock core of mass 0.912M E 
surrounded by a frozen, salty, H 2 0 ice mantle of thickness 150 km seems to do a very good 
job at explaining the Galileo gravity data. 


Table 9: Structural Properties of the Computed Satellite Models 


Model 

Europa 

Ganymede 

Ganymede 

Callisto 

Callisto 

Property "" — 

[2-zone] 

[2-zone] 

[3-zone] 

[1-zone] 

[2-zone] 

Physical radius (km) 

1565 

2634 

2634 

2403 

2403 

Central pressure (kbar) 

40.7 

74.4 

94.0 

31.6 

34.5 

Central temperature (K) 

1073 

1073 

1073 

253 

263 

Central density (g/cm 3 ) 

3.776 

3.888 

5.607 

2.142 

2.326 

Rock core mean 

680 

820 

820 

- 

200 

temperature (K) 





154 

Core edge temperature (K) 

151 

209 

209 

- 

Core radius (km) 

1415 

1640 

890/1635 

- 

2137 

Convective boundary (km) 

1492 

2172 

2174 

2285 

2230 

Mean density (g/cm 3 ) 

2.989 

1.888 

1.896 

1.860 

1.803 

Moment-of-inertia factor 

0.3503 

0.3156 

0.3106 

0.3836 

0.3580 
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EUROPA MODEL 


DEHYDRATED 
ROCK CORE 

• 1415 km radius 

• Mass: 0.9116 M E 

• Density 

(1 -bar, 298.15 K): 
3.68 g/cm 3 



FROZEN, SALTY 
H 2 0 MANTLE 

• Thickness: 150 km 

• Mass: 0.0884 M E 


• Ice 1 h/lce 2 
boundary: 
1457 km 


1565 km 


Model mean density: 2.99 g/cm3 

Predicted moment-of-inertia coefficient: 0.350 ± 0.002 

Figure 5. Cross-sectional view of the proposed 2-zone model for Europa. A frozen mantle of salty H 2 0 
ice of thickness 150 km and mean temperature 140 K surrounds a dehydrated rock core of 
mean temperature 680 K. NaCl is the principal salt constituent of the ice. The chief 
constituents of the rock are FeS, Fe 3 0 4 , MgSiOg and Mg 2 Si0 4 . 


8.3 Ganymede 

Consider now Ganymede. Gravity observations by the Galileo spacecraft have 
confirmed the long-held suspicion that this satellite is a heavily differentiated structure, 
consisting of a dense rock/metal core surrounded by a thick mantle of H 2 0 (Anderson et 
al., 1996a). The gravity data can be accounted by a 2-zone model provided that the 
uncompressed rock has a density of order 3.6 -3.7 g/cm 3 . But the unanticipated discovery 
of an intrinsic Ganymedan magnetic field by the Galileo magnetometer experiment 
(Kivelson et al., 1996) favors a 3-zone model in which Fe and FeS have separated out 
from the rock to form an inner metallic core. If such a core were molten, then dynamo 
action within it may generate the observed magnetic field (Schubert et al., 1996). 

Prentice (1996c) noted that 2-zone models of Ganymede also have the capacity to 
produce a native magnetic field if one takes into account a process of thermoremanence 
within the magnetite constituent of the rock. That is, the release of radiogenic heat during 
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the course of Ganymede’s evolution causes the temperature in certain parts of the rock core 
both to rise above and then descend again through the Fe 3 0 4 Curie point [T c = 858 K]. If 
this event happened early in the satellite’s evolution, at a time when the ancient Jovian 
magnetic field was, say, 10 times larger than it is today, and as long as this field also 
remained stable in direction over long periods of time —10^ yr, then thermoremanent 
magnetization (TRM) of the cooling magnetite may have occurred. The reader is referred 
to Dunlop and Ozdemir (1997, p. 5), for an excellent description of TRM. Quite 
independently, Crary and Bagenal (1996) have also proposed that remnant magnetization 
may be responsible for Ganymede’s magnetic field. In a more recent paper (Crary and 
Bagenal, 1998), these authors favor a paleomagnetic field arising from dynamo action 
within a previously molten core rather than a background Jovian field as the source of 
magnetization in Ganymede. 

In the present paper it is assumed that Ganymede’s observed magnetic field arose 
from TRM fed by an early, powerful, magnetic field of Jupiter. The calculations are there- 
fore confined to 2-zone models. Further, it will be assumed that the initial Ganymede (G) 
had already undergone complete differentiation between its rock and ice phases at the time 
of its formation, through the disposition of substantial amounts of gravitational accretion 
energy (Lunine and Stevenson, 1982). Next, as for Europa, it is assumed that the release of 
radiogenic heat causes a complete dehydration of the rock core, including the eviction of 
NaCl. From Table 4 and 6, it follows that the mass of the rock core is 0.4786 Mq, while 
the mass of the ice mantle includes 0.04555 Mq water of hydration and 0.00075 Mq of 
NaCl, in addition to the native ice mass of 0.4751 Mq. That is, the Ganymedan ice is only 
9% as salty as that of the Europa’ s ice mantle. Returning to Table 4, the new mass 
fractions of MgSi0 3 , Mg 2 Si0 4 , MgAl 2 0 4 , Si0 2 , and Na-K feldspars, relative to the initial 
rock mass, are 0.2325, 0.1759, 0.0097, 0.0000 and 0.0689, respectively. The RTP density 
of the dehydrated salt-free rock core is 3.6758 g/cm J . 

Figure 6 gives the temperature distribution T(r,t ) within Ganymede versus internal 
radius r at various times t in the course of its evolution to present age 4600 Myr. After 
-150 Myr. the temperature at the rock/ice boundary at radius r x = 1640 km rises to the 
value -209 K where the ice becomes convective at the pressure p\ = 20.2 kbar at that point. 
The thickness of the convective layer of ice then grows progressively to -530 km by 
4600 Myr. Unlike the case for Europa, the normal creep criterion given by equation (9) 
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does apply to Ganymede since dT me i t (p)/*7p>0 in the relevant pressure regime. The 
pressure and temperature at the outer edge of the convective layer are 7.95 kbar and 172 K. 

Most of the main parameters of the computed 2-zone structural model for Ganymede 
are listed in Table 9. The outer mantle of H 2 O consists of an inner convective region of 
mass 0.2453 M g and mean temperature 190 K and an outer conductive zone of mass 
0.2760 M g and mean temperature 130.5 K. The H 2 0 ice VIII/VI phase boundary lies 
within the convecting region at the 1715 km radius point. The VI/II and Il/Ih phase 
boundaries are at radii 2200 km and 2556 km. 



Figure 6. Internal temperature profiles of a 2-zone model for Ganymede versus fractional radius 
x = r/r G plotted at several times during the evolution to age 4600 Myr. The surface 
temperature of the satellite is assumed to be constant (i.e., 103 K) throughout the 
evolution. After time ~ 150 Myr, the temperature T, at the grid-point radius x l = 0.635 of 
the rock-core/ice-mantle interface rises to 0.6 of the local H 2 O ice melting temperature 
^melt = 345 K - Thereafter, the ice becomes convective and the temperature at the rock 
core boundary stays fixed. 
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Comparing the mean density p of the 2-zone satellite model with the observed value 
Pq =1.936 ± 0.022 g/cm 3 (see Table 5) we see that p falls short of p G b y 0 048 g/cm 3 , 
or 2.2 o. Next consider the mean axial moment-of-inertia factor/. Anderson et al. (1996a) 
report a value / G = 0.3 105 ± 0.0028. Hence the 2-zone model value of 0.3 156 exceeds f G 
by ~1.8 a. Thus although the parameters of the 2-zone model are not in serious 
disagreement with the Galileo data, neither are they satisfactory. Suppose now we relax 
the 2-zone model assumption and introduce a dense inner core of solid FeS at the center of 
the satellite. This can be achieved by melting the FeS constituent of the rock (T me j t = 

1468 K) early in the satellite’s radiogenic evolution. A preliminary structure calculation, 
based on the same thermal profile as the 2-zone model, indicates that if all of the FeS is 
brought to the center, then the axial moment-of-inertia factor falls to 0.3106. This is 
almost the same as/ G . The density of this 2-zone model increases slightly to 1 .896 g/cm 3 
but it still falls short of the observed value by -0.040 g/cm 3 . 

To achieve a model for Ganymede that is consistent with both the observed mean 
density and axial moment-of-inertia, it is necessary to increase the rock mass fraction of the 
satellite. According to the condensation model for the Galilean satellites presented in 
Section 4, the rock mass fraction of Ganymede is controlled by the enhancement factor Cen- 
This is the abundance of heavy elements in the primitive Jovian atmosphere, relative to 
solar abundances. In this paper, the value Q n = 2.08 was chosen to match the H 2 S to H 2 
mixing ratio of 7.7 x 10‘ 5 observed in the Jovian atmosphere by the Galileo Probe 
(Niemann et al., 1996, 1998). If the H 2 S to H 2 ratio in the deep atmosphere is -9 x 10 5 , 
as suggested by the data in Figure 6 of Niemann et al. (1998) so that £ en - 2.5, the mean 
density of Ganymede can be accounted by the condensation model presented here. The 
resulting 3-zone model for Ganymede would consist of an inner dense core of solid FeS of 
mass -0. 1 16M g , surrounded by an FeS-free rock layer of mass 0.382 Mq and an outer, 
frozen, salty ice mantle of mass ~0.502M G . This model will not be pursued any further 
here. The conclusion of Anderson et al. (1996a) that Ganymede is a 3-zone structure 
containing an inner metallic core is thus a robust one. 

Lastly, we return to the origin of the Ganymedan magnetic field. Figure 7 shows a 
plot of the temperature at the 1410 km and 1475 km radius points as a function of time. 
During the radiogenic evolution to present age, the temperature of the rock between these 
points both rises above and descends through the FejC^ Curie point. As discussed above, 
the magnetite constituent of this 65 km thick layer may thus become permanently 
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magnetized through thermoremanent magnetization fed by Jupiter’s ancient magnetic field. 
Most of the field acquisition took place at a time ~2 x 10 9 yr. Jupiter’s magnetic field must 
have then been at least one order of magnitude larger than the present day ambient field if 
we are to account for the strength of the Ganymedan dipole (Kivelson et al., 1996). Figure 
8 is a simplified schematic picture of Ganymede, showing the magnetized layer and other 
principal features of the 2-zone model. 



TIME f/10* YR 


Figure 7. The temperature at the fractional radius grid points x = r/R c = 0.535 and r/R Q = 0.56 within a 
2-zone model of Ganymede plotted as a function of time. Between these two radius points, 
which are located in the dehydrated rock core, the temperature rises above and then descends 
through the Fe 3 0 4 Curie point during the thermal evolution to present age. Ganymede’s 
observed magnetic field may thus have been created by thermoremanent magnetization fed by 
Jupiter’s ancient, but powerful, magnetic field (Prentice, 1996c). 
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A POST-GALILEO, TWO-ZONE 
GANYMEDE MODEL 


DEHYDRATED ROCK CORE 

• Mass: 0.479 Mq 

• Radius: 1640 km 

• Mean Density: 3.74 g/cm 3 

• Mean Temperature: 820 K 

• Central Temperature: 1073 K 

• Central Pressure: 74.1 kbar 


MODEL CHARACTERISTICS 

• Mean Density: 1.89 g/cm 3 

• Moment-oMnertia 
Factor: 0.316 ±0.002 



PURE WATER ICE MANTLE 

• Mass: 0.521 M G 

• Thickness: 994 km 

• Mean Density: 1 .30 g/cm 3 

• Mean Temperature: 160 km 
Surface Temperature: 103 K 


MAGNETIZED Fe 3 0 4 
SHELL: 65 km thick 


2634 km 


Figure 8. An illustrative summary of the 2-zone structural model for Ganymede. A frozen mantle 
of H 2 0 ice of mass 52.1% of the Ganymede mass surrounds a dehydrated rock core. 
Dissolved NaCl makes up 0.14% of the ice mass. Ganymede’s native dipolar magnetic 
field arises from a 65 km thick shell of rock containing magnetized magnetite. 


8,4 Callisto 

It remains to consider Callisto (C). According to Figure 3 and Table 4, the 
condensate from which the initial Callisto formed consisted of an intimate mixture of 
hydrated rock, H2O ice and NH3 H2O monohydrate. In a pre-Galileo analysis, Prentice 
(1996a) observed that for the choice of heavy element enhancement factor ^ en = 2, the 
masses of rock and ice in Callisto occur in just the right proportions to produce a 
chemically uniform satellite model whose mean density p matched the observed value 
p c =1.852 ± 0.012 g/cm 3 (see Table 5). On the basis of this good agreement, it was 
predicted that the satellite had remained cold, structurally unmodified and magnetically 
inert since the time of its formation (Prentice 1996c). Figure 9 provides a schematic view 
of this pre-Galileo model of Callisto. The basic structural parameters of the model are also 
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PRE-GALILEO 
CALLISTO MODEL 


CHEMICAL COMPOSITION 
UNIFORM MIXTURE OF: 

• Hydrated Rock (49.6%) 

• Water Ice (45.7%) 

• Ammonia Ice (4.7%) 
(Present as NH 3 *H 2 °) 


ICEIh 


MODEL CHARACTERISTICS: 

• Mean Density: 1.66 g/cm 3 

• Central Temperature: 253 K 

• Central Pressure: 31.6 kbar 

• Surface Temperature: 103 K 

• Moment-of-lnertia factor: 0.364 ± 0.002 



ICE CONVECTION 
BOUNDARY 
2285 km 


2403 km 


ICE 1h/2 BOUNDARY 
2338 km 


PREDICTED 

MAGNETIC STRUCTURE: 

MAGNETICALLY 

INERT 


Figure 9. Cross-sectional view of the internal structure of a pre-Galileo Callisto model consisting of a 
uniform mixture of hydrated rock, H 2 0 ice and NH 3 ice. The dashed circles mark the 
boundaries between the various high pressure phases of H 2 0 ice. The model remains cool 
and hence magnetically inert throughout its entire radiogenic evolution through efficient 
solid-state creep of the ice. 


listed in Table 9. In specifying an equation of state for NH 3 H 2 0 it is assumed, for 
simplicity, that this substance behaves like a simple molar mixture of H 2 0 and NH 3 ices. 
According to Croft et al. (1988), the density of NH 3 ice at 1-bar pressure varies with tem- 
perature according as Pnh 3 = 0.8659<T' ?(7) g/cm 3 , where g(T) = 4.198 x 10~ 7 r 2 - 2207 . 
The compressibility of NH 3 ice was modelled with data from Stewart (1960). 

It was reassuring to learn that the initial gravity measurements of Callisto obtained 
during the C3 flypast on November 4, 1996, indicated that this satellite does appear to have 
a chemically undifferentiated structure (Anderson et al, 1997b). Furthermore, the Galileo 
magnetometer observations indicated the absence of any internal magnetic field (Khurana 
et al, 1997). More accurate gravity data gathered by the Galileo Orbiter during the CIO 
flypast of September 17, 1997, has however, produced a substantially revised picture for 
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the internal structure of Callisto (Anderson et al., 1998b). Instead of being chemically 
homogeneous, as the earlier data indicated, it now seems certain that a partial separation 
between the rock and ice phases has occurred. The value for the mean axial moment-of- 
inertia factor f c obtained from the CIO gravity data is 0.358 ± 0.004. This lies well below 
the value -0.384 expected for a chemically uniform rock/ice satellite (see Table 7). Yet 
despite the gravitational evidence for differentiation, the magnetometer observations still 
indicate the absence of any central, native magnetic field (Khurana et al., 1998). Most 
likely, therefore, the core of the satellite is still cold and thermally inactive, as predicted 
earlier Prentice (1996a, c), and ice is mixed with the rock all the way to the center. 

A new post-Galileo Callisto model has been constructed in which it is assumed that 
gravitational accretion energy, which had previously been ignored, has melted all the ice in 
the outer 34.5% of the mass of the primitive Callisto. The rock component of this melted 
outer layer then migrated into the frozen inner core. The resulting satellite structure 
consists initially of a rock-rich, icy core of total mass 0.826 Mq , surmounted by a liquid 
mantle of H 2 0 and NH 3 . The liquid froze fairly quickly to form an icy mantle of thickness 
-270 km. The mass fractions of hydrated rock (the same mixture as in Table 6), H 2 0 ice 
and NH 3 ice of the core 0.6000, 0.3623 and 0.0377, respectively. The 1-bar density of the 
core material at the typical temperature of 200 K is 1.593 g/cm 3 . The rock volume fraction 
^rock = 0-3° 5 of the core material is less than the critical value -0.4 at which the rock/ice 
mixture becomes too viscous to efficiently convect (Friedson and Stevenson, 1983). We 
can safely assume therefore that: (i) solid state convection will continue to operate 
throughout the satellite during the course of its radiogenic evolution to age 4600 Myr, and 
(ii) that this convection renders the core chemically uniform. 

The structural properties of the new, partly differentiated, model for Callisto are 
given in Table 9. Figure 10 provides a simple schematic of this satellite model. We 
observe with some excitement that the axial moment-of-inertia factor of this first-try 
model, namely /= 0.358, coincides exactly with the value/ c derived from the Galileo CIO 
flypast (Anderson et al, 1998b). Clearly the choice X rock = 0.60 for the core rock mass 
fraction proved to be the right one. A larger X rock would have resulted in a smaller/. The 
central density of the satellite, viz. 2.33 g/cm 3 , agrees well with the mean core density of 
2.3 g/cm 3 estimated by Anderson et al (1998b) for their two-layer model. This excellent 
agreement provides a valuable independent check on the validity of the computational 
scheme for calculating satellite stmctures which has been presented here. 
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ROCK-RICH ICE CORE 

MASS: 0.826 M c 

CHEMICAL COMPOSITION: 

• Hydrated Rock (60%) 

• Water Ice (36.2%) 

• Ammonia Ice (3.8%) 


A POST-GALILEO 
CALLISTO MODEL 


CORE/MANTLE 
BOUNDARY 
2137 km 



PURE ICE MANTLE 

MASS: 0.174 M c 

CHEMICAL COMPOSITION: 

• Water Ice (90.6%), 

• Ammonia Ice (9.4%) 


MODEL CHARACTERISTICS: 

• Mean Density: 1.80 g/cm 3 

• Central Temperature: 263 K 

• Central Pressure: 34.5 kbar 

• Central Density: 2.33 g/cm 3 

• Moment-of-inertia factor: 0.358 ± 0.002 


2403 km 


H 2 0 ICE 1 h/2 
BOUNDARY 
2291 km 


PREDICTED 

MAGNETIC STRUCTURE: 
MAGNETICALLY 
INERT 

CONFIRMED 


Figure 10. Schematic view of the cross-section of a post-Galileo, partly differentiated model for 
Callisto’s internal structure. A mantle of pure H 2 0 and NH 3 ices of mass 0.174 times the 
satellite mass surrounds a denser core containing a uniform mixture of hydrated rock and 
ices. Solid-state convection of the ice keeps the satellite cool to the center. The first 90 km 
of the outer ice mantle is also convective. Although the moment-of-inertia factor of this 
model matches the Galileo observation (Anderson ct ul. , 1998b), the mean density of the 
model falls short by 0.05 g/cm 3 . A further improvement in modelling is required, especially 
to allow for NHyHiO formation, as discussed in Section 9. 


There is, unfortunately, one major deficiency with one aspect of the proposed new 
Callisto model, namely that the calculated mean density p is 0.049 g/cm 3 less than the 
observed value, namely p c =1.852±0.012 g/cm 3 . This happens because the ice is now 
less compressed than that of the previous homogenous model, since not all of it is subject 
to the weight of overlying rock. Yet quite remarkably, this shortfall in the p of the new 
Callisto model is exactly the same amount as the density shortfall in the Ganymede model 
presented in Section 8.3. Both satellite models could, therefore, be corrected 
simultaneously by choosing a larger value for the enhancement factor C en of heavy 
elements in the primitive Jovian atmosphere, as was discussed in that Section. There is 
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probably no point, however, in carrying out such a refinement in the satellite modelling 
until all the Galileo data is in. In this regard, it is appropriate now to take into account 
improved determinations of the satellite mean radii (Davies et al., 1998; Thomas et al., 
1998). The calculations reported here have been based on the Voyager determinations 
(Davies et al., 1992). 


9. Conclusions and Discussion 

We have considered a reformulated version of the original nebula hypothesis of 

Laplace which includes the influence of turbulent stress (^p,v, ^ = f3pGM{r)l r arising 
from supersonic convective motions within each of the primitive contracting clouds which 
formed Jupiter and the Sun. This ‘modem Laplacian theory’ of solar system origin has the 
capacity to account simultaneously for the observed broad physical and chemical structure 
of both the planetary and Galilean satellite systems for a single choice of the turbulence 
parameter (3, viz. [3 ~ 0.085. Furthermore, several key pre-Galileo predictions for the 
properties of the Jupiter system which were made on the basis of this theory have now been 
confirmed by the Galileo spacecraft data. These predictions include a 2-fold enhancement 
of the sulphur abundance in the Jovian atmosphere relative to the solar level of that 
element, a mean axial moment-of-inertia factor for Europa which coincides exactly with 
the Galileo Orbiter measurement, and the absence of any native magnetic field of Callisto. 

We have suggested that each of the Galilean satellites condensed individually from a 
concentric family of orbiting gas rings. These rings, each of about the same mass 
~ 1.6 x 1 0 25 kg, were cast off by the contracting proto-Jovian cloud close to the present 
satellite orbits. Setting aside the heat due to the primitive sun, the temperatures T n of the 
gas rings vary with mean orbital distance R n according as T n = 197(/?q !R n ) K , where 

R g denotes Ganymede’s present distance. This law, along with a similar scaling law for 
the gas pressures p n on the mean orbits, yields a suite of bulk chemical compositions of the 
condensate in each gas ring whose mean densities closely mimic the observed trend of 
mean densities of the Galilean satellites. Although Europa condensed above the H 2 0 ice- 
point, the condensate rock at its orbit contains a large quantity of the highly hydrous 
mineral brucite. Early decomposition of this mineral through radiogenic heating created 
Europa’s present -150 km thick ice mantle. The condensate at Callisto’s orbit includes 
4.75% by mass of NH 3 . 
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Detailed thermal evolutionary and physical structural models for each of the ‘icy’ 
Galilean satellites were constructed on the basis of the suite of bulk chemical compositions 
which follow from the above condensation model. By constructing both homogenous and 
chemically differentiated models, it is possible to infer which satellite structure best 
explains the Galileo spacecraft observations. A 2-zone model for Europa, consisting of a 
dehydrated rock core of mass 0.912 M E and a 150 km thick outer mantle of frozen, salty, 
H 2 O ice is consistent with all of the gravity data. For Ganymede, a similar 2-zone model, 
based on the initial 52.5% : 47.5% mass proportions of hydrated rock and ice which follow 
from the condensate model for the heavy-element enhancement factor £ en = 2.08 in the 
proto-Jovian cloud, was constructed. This yields a structure whose mean density and axial 
moment-of-inertia factor are in error by ~2 a from the observed values. A 3 -zone model of 
Ganymede, containing an inner core of solid FeS and a greater total mass fraction of rock, 
is needed to match all of the data. This may be achieved by assuming a higher value for 
the heavy element enhancement factor <^ en , namely £ en « 2.5. 

It is suggested that Ganymede’s native dipole magnetic field was created by 
thermoremanent magnetization some ~2 x 10 9 yr after this satellite was formed. This 
activity took place within a cooling shell of magnetite-bearing rock of radius ~1450 km 
and thickness ~65 km. Thermoremanence was fed primarily by Jupiter’s early, powerful 
magnetic field (Prentice 1996c). Ancient dynamo action within the inner core of FeS, 
which was initially in a molten state, may also have assisted in the field acquisition 
process. 

The present work suggests that Callisto is a cold, thermally dormant, and 
magnetically neutral structure with H 2 0 and NH 3 ices extending all the way to the satellite 
center. The Galileo Orbiter measurements indicate, however, that Callisto is more 
centrally condensed than a chemically homogeneous body. The gravity data can be 
accounted for by a modified 2-zone model consisting of a rock-enriched icy core of mass 
0.826 M c (60% hydrated rock, 40% ice) surrounded by a pure H 2 0-NH 3 ice mantle of 
thickness ~270 km. Such a mantle may have formed through melting of the outer satellite 
layers in the late stages of accretion. If the same process also occurred with Ganymede, the 
larger mass of this satellite coupled with its higher mass fraction of rock would result in 
more extensive melting and hence a higher fraction or rock precipitating into the core. 

Such an event would serve to stymie efficient solid-state convection during the satellite’s 
subsequent radiogenic evolution and so lead to a completely differentiated structure. Thus 
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it was a combination of accretional melting (Lunine and Stevenson, 1982) coupled with a 
subsequent loss of convective efficiency (Friedson and Stevenson, 1983) which, most 
likely, determined the different pathways followed by Ganymede and Callisto. In short, 
Ganymede was ‘too big for its boots.’ Callisto’s capacity to stay cool was also assisted by 
this moon’s greater mass fraction of ice. 

Lastly, I note that the very recent discovery by the Galileo magnetometer experiment 
(Khurana et al., 1 998) of a variability in the Jovian magnetic field near Callisto may be 
accounted for by the proposed store of NH 3 . The refreezing of Callisto s initial liquid 
outer mantle of H 2 0-NH 3 results in the formation of a eutectic mixture of H 2 0 and 
NH 3 -2H 2 0 ices (J.S. Kargel, priv. comm). But the melting temperature of NH 3 -2H 2 0 ice 
is very low, namely T m (0) = 176 K at pressure p — 0 kbar (Hogenboom et al., 1997). It is 
possible, therefore, that radiogenic heat arising from Callisto’s rock-ice core may be just 
sufficient to maintain a narrow shell of liquid NH 3 -2H 2 0 at the pressure depth p -3 kbar. 
This depth, where T m achieves a local minimum -177 K, lies just outside the dense rock- 
ice core. Dissolved NaCl in this liquid ammonia dihydrate ocean would then create the 
electrically conducting layer inferred by the Galileo observations. It was Lewis (1971) 
who first pointed out that NH 3 may lead to a liquid ocean in Callisto. The present work 
suggests that Callisto alone possesses an internal ocean. Europa and Ganymede are 
expected to be frozen solid and Europa’s magnetic field, if one exists (Kivelson et al., 
1997), arises through thermoremanent magnetization (Prentice, 1996c, 1999). 

The Galileo Project has revealed that the Galilean satellite system of Jupiter is 
remarkably diverse in its physical and chemical structure. The observations have provided 
valuable clues to the conditions under which this system formed. The model of satellite 
origin presented here has the capacity to quantify these conditions and appears to go a long 
way towards explaining much of the data. Hopefully this work will make a lasting 
contribution towards understanding how our solar system came into being. 
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